function  interp_function  = Interp_function(estim_value,log_Skills_tp1,H_prime)

interp_function_hat(:,:,:,1)= estim_value(1).*ones(size(log_Skills_tp1));
interp_function_hat(:,:,:,2)= estim_value(2).*log_Skills_tp1;
interp_function_hat(:,:,:,3)= estim_value(3).*H_prime;
interp_function_hat(:,:,:,4)= estim_value(4).*log_Skills_tp1.*H_prime;
                             
interp_function = sum(interp_function_hat,4);          